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The simple model of a semiconductor, based on a single effective mass for 
the charge carriers and a spherical shape for the surfaces of constant energy, 
is now Jcnoivn to he inadequate for most of the semiconductors which have 
been extensively studied experimentally. However, some of these do corre- 
spond to what may be called the "many-valley" model, a model for which 
the band edge occurs at a number of equivalent points K*'^ in ivavc number 
space, and for which the surfaces of constant energy are multiple ellipsoids, 
one centered on each of these points. This paper develops, for models of this 
type, the theory for: mohility (Section 2) and its temperature dependence 
{Section 3); thermoelectric power (Section 4); piezoresistance (Section 5); 
Hall effect (Sections 6 and 9); high-frequency dielectric constant (Section 7); 
and magnetoresistance (Sections S and 9). These phenomena are treated, 
for cases to which Maxwcllian statistics apply, on the assumption that the 
scattering of the charge carriers is describable by a relaxation time which de- 
pends on energy onhj, but is otherwise unrestricted. This assumption can 
he shoivn to he justified in a large class of cases, although for some cases it 
fails, notably when ionized impurity scattering predominates and at the 
same time the effective mass is very anisotropic. Special attention is given to 
the role of inter-valley lattice scattering, i.e., to processes whereby a charge 
carrier is scattered from the neighborhood of one of the band edge points 
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K^'' to the neighborhood of a different one. Numerical calculations are pre- 
sented which show the effects of such processes on the magnitudes and tem- 
perature variations of the effects listed above. 

1. THE MANY-VALLEY MODEL 

Most of the literature of semiconductor theory has been based on what 
we shall call the simple model. This model is based on the assumption 
that the minimum energy in the conduction liand, or the maximum 
energy in the valence band, is possessed by only ont^ (juantum state of 
either spin. This state has the form of a Bloch wave with wave number 
K = 0.* States with energies near the band edge value therefore have 
small K values, and, since their energies t{K) must vary continuously 
with K in this region, e(K) for small K must be a quadratic form in 
Ki , Ky , K^ Ai the crystal stiiicture is cubic, e(K) a K\ and the sur- 
faces of constant energy are spheres in K-space. 

It has long been known that other models are possible, and indeed 
likely in many cases. In recent years it has become clear that the simple 
model does not apply to any of the four cases corresponding to n- and p- 
type germanium and silicon. The evidence for this includes magneto- 
resistance^' ^' ^ and piezoresistance^ effects, cyclotron resonances,^ and 
many other phenomena. Now the possible alternatives to the simple 
model are the various models for which there is more than one state, 
apart from spin degeneracy, with the band edge energy. These models 
fall into two general categories. 

(A) Models for which the band edge energy occurs for several wave 
number vectors K'*', but for which there is only one state of each spin 
having this energy and a given K''^. For a conduction band model of 
this sort the energy e, considered as a fimction of K, has a number of 
minima or "valleys", hence we shall call these "many-valley" or "simple 



* For the convenience of the reader the notations defined in the text are re- 
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many -valley" models. For a valence band the situation is similar, but 
in\'erlod. 

(B) Models for which, apart from spin degeneracy, there are two or 
more states with the hand edge energy and the same wa\'e number 
vector. These we shall call "degenerate" models. We may subdivide them 
further into "degenerate single-valley" and "degenerate many-valley" 
cases according to whether the band edge energy occurs for only one 
wa\'e vector, or for several. 

Tliis paper is concerned with the transport properties of the simple 
many-vall(\v models defined under (A). These models are much simpler 
to handle than the degenerate types, for reasons which are illustrated 
in Fig. 1. This illustration shows schematically the form in wave number 
space of the surfaces of constant energy, near the band edge energy, for 
four models. For the simple model, shown in (a), the locus of a given 
energy is, as already stated, a sphere. For a simple many-valley model, 
shown in (b), the locus is a set of ellipsoids centered about the band 
edge points K'*'. The ellipsoidal shape is required by the facts that 
energy must depend contiiiuously and different iably on K and have an 
extremum at each K*'\ For a degenerate model, however, the dependence 
of energy on the components of K is singular at the band edge point,' 
in that unique second derivatives do not exist: energy varies quadrati- 
cally wiMi K in any given direction from this point, but the coefficients 
going with different directions are determined by a secular equation. 
The result is that the contours of constant energy may look as shown in 
Fig. 1(c) (degenerate single-valley case). Degenerate many-valley cases 
are of coui'sc similar, but with the surfaces multiplied, as in Fig. 1(d). 
Such situations are obviously harder to handle mathematically than 
those of Fig. 1(b). 

Hesirlcs the irregularity of the energy surfaces, there is another dif- 
ference between these two types of cases which greatly complicates theo- 
retical work with degcneiate models. This is that in most cases the 
energies of the two or more states going with a gi\'en band-edge K will 
be split by spin ori)it coupling. If this splitting is «A:'y' it can usually be 
ignored, and if it is »/.■?' it may elTectively convert the degenerate model 
into a simple or simple many-valley model, rnfortunately, it usually 
happens that neither of these extremes applies, and for such inter- 
mediate cases not only do we have to deal with energy surfaces like 
those of Fig. 1 (a), l)ut, much worse, the variation of energy with K is not 
a simple ([Uadratic dependence even in a fixed directi{)n from the band 
edge point. 

In view of all these complications of the degenerate models, it is for- 
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tunate that the shnple many-valley case does seem to occur for n-type 
silicon and n-type germanium. One of the best ways of finding out 
whether it occurs for any given semiconductor is to compare obsei-va- 
tions on this material with theoretical predictions for the various possible 
models of the simple many-valley type. We proceed now to derive these 
predictions, assuming for simplicity that the charge carriers have 
Maxwellian statistics and an effective relaxation time dependent only 
on energy. We shall take up the simplest properties first, the more com- 
plicated ones later. Sections 2 to 5 will consider perturbation of the 
distribution function of the carriers by a static electric field, Section 6 
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Fig. 1 — Different types of band structure for a semiconductor, illustrated by 
the forms of the surfaces of constant energ>' in wave number space. The band 
edge points are represented by heavy dots. 
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the Hall effect, Section 7 the perturbation by an oscillating electric 
field, Section 8 magnetoresistance, and Section 9 effects at high mag- 
netic fields. 

In presenting this material our primary objective will be to provide a 
coherent treatment of all the effects in language as simple and physical 
as possible. Thus, for example, the Hall and magnetoresistance effects 
will be discussed ab initio, although many of the details presented here 
have been derived and published independently by several workers.^' ^ 
Nor is the theory of this paper the ultimate in refinement: at cost of a 
little more mathematical complication, the assumption of a relaxation 
time dependent only on energy can be dispensed with, and anisotropy in 
the scattering processes acting on the carriers can be taken into ac- 
count.'' However, the present simpler treatment illustrates most of the 
physical principles involved in the various phenomena, and turns out to 
be quantitatively adequate in a large class of cases. 

2. CONDUCTIVITY 

In this section we shall solve the Boltzmann equation for the effect of 
a constant electric field E on the motion of charge carriers in a simple 
many- valley band. Maxwellian statistics will be assumed. Thus if 
AP = /t (K — K^'') measures the deviation in crystal momentum 
space from one of the band edge points K''^, then f or E = the proba- 
bility of occupation of the state described by AP (by an electron or hole, 
depending on the sign of the carriers) is 



■(0) 



r = exp 



2mi* 2m-* 2m3* 

It 



(1) 



where tf is the Fermi level, ei, the band edge energy and mi*, fn2*, ms* 
are the effective masses in the three coordinate directions 1, 2, 3 which 
are principal axes for the energy surfaces of the valley in question. 
When E = F^ the distribution function / is determined by the competi- 
tion between the perturbing effect of E and the restoring effect of scatter- 
ing processes which try to restore the form (1). To make the problem 
tractable we shall assume that the scattering processes which the charge 
carriers undergo are described by a relaxation time which is a function 
of energy e only. In other words, we shall assume that for any slight 



* C. Herring and K. Vogt, to be published. 
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departure of /from/'*" the time rate of change of /due to collisions is 

df\ _ (/ - r) (2j 



dtj. Tie) 

The legitimacy of this assumption is analyzed in Appendix A. It is shown 
there that the assumption should be rather accurately valid for all kinds 
of inter-vatley scattering — defined as scattering from the neighborhood 
of one band edge point K*'' to the neighborhood of another K^^' — and 
for intra-valley lattice scattering due to optical modes or to neutral 
impurities, provided, in the latter case, that the temperature is low 
enough. For intra-valley lattice scattering due to acoustical modes the 
assumption t = T(e) is not necessarily valid, but the arguments of Ap- 
pendix A suggest tht it will often be a good approximation. For scatter- 
ing by ionized impurities, however, this assumption will usually be a 
poor approximation. There is a good prospect that in the near future 
the adequacy of this approximation for lattice scattering can be quanti- 
tatively estimated for some substances. If it should turn out to be in- 
adequate, the necessary generalization of the calculations of this paper 
can probably be made without great effort. 

With the assumptions just stated in (1) and (2), the Boltzmann equa- 
tion for a steady state in the presence of a field E takes the form 

0=%- ±.E-V./- ^■^"■^""^ (3) 

dt T 

where the upper sign is for electrons in a conduction band, the lower 
for holes in a filled band. If, as is customary, we set 

/=.r +E-f'" -\-0{E'), (4) 

(3) gives, just as in the simple theory, 

f"' = ±erVrf' (5) 

Having obtained the solution of the Boltzmann equation in the foi-m 
(4), (5), we shall now evaluate the electron current density j fi'om it. If 
/"*' is Maxwellian, 

V,f<«> = 'V^ V.A. = -^ r (6) 

c/Ae kl 

where v is the group velocity and Ae ^ | e — fb | is the distance from the 
band edge. The contribution of carriers in the ?th valley to the current 
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density is then 



j"' - E (±e)v(AP*''')/(AP"') 

AP('),s 

- 7^ E rr{Ae)E-w (7) 

where the summations are over all AP''' occurring in the iih valley in 
unit volume of material, and over both states of spin. 

The expression (7) states that any single valley i possesses an aniso- 
tropic conductivity tensor o-o^''', or an anisotropic mobility tensor jUa^^'*, 
i.e., 



where 



n 



= Er (9) 



ip(n,. 



is the number of carriers in valley i per unit volume. If we choose the 
.r, //, and z axes to be along the principal axes of the ellipsoidal energy 
surfaces of valley /, (7) shows that o-„^ ' and ^03 '* will be diagonal. Each 
diagonal element ^u''' will involve a Maxwellian average of v>,^T(Ae). 
Now the equipartition principle leads us to expect that the average, over 
an energy shell Ae to Ae + c/Ae in AP-space, of the kinetic energy asso- 
ciated with the a;-component of velocity should be the same as that as- 
sociated with the y- or s-component. This is easily demonstrated ex- 
plicitly (Appendix B). Thus, if mi*, nh*, nis* are the effective masses 
in the three principal directions, 

H nn*v{ - H m.^v:' = ^ m,*v,^ (10) 

Tiierefore (7) and (8) give, in our system of axes, 



(i) 



e <Aer> 



ilea, ^ ^ 

(i) 



M.8 =0 (a ?^ /3) (12) 

where the angular brackets denote Maxwellian averages: 

<Aer> - ll^^rr/T.r\^t^. (13) 

^p ip 

The denominator of (11), <Ae>, of course equals % kT by equipartition. 
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The formula (11), it will be noted, is the same as that of the simple 
theory^ with m* replaced by ??;„*. 

The overall conductivity tensor due to the carriers in all the valleys is 
of course SiO-a,?''', and the overall mobility tensor is the average of 
Mop''' over the different valleys. For a cubic crystal the mobility is the 
same in all directions, so we have 



ZNy . . (14) 

„ e <A€T> 
~ my> <Ac> 

where Ny is the number of valleys and m' is an average inertial mass 
defined by 

1 






\_mi* m-t m3*. 



(15) 



This mass, as we shall see in Section 7, is the one most directly meas- 
ured by the Benedict-Shockley experiment on high-frequency dielectric 
constant. 

3. TEMPEEATURE VARIATION OF LATTICE MOBILITY 

The T occurring in the mobihty expression (14) differs from the r of 
the simple model in that it contains the effect of inter -valley scattering 
in addition to intra-valley and impurity scattering. Inter- and intra- 
valley scattering differ in that most of the phonons emitted or absorbed 
in intra-valley scattering have energies « the energies of the charge 
carriers, while those involved in inter-valley scattering usually do not. 
If K''' and K'''' are two different band edge points, scattering of a carrier 
from valley i to valley y must involve emission or absorption of a phonon 
of wave number close to iq.y , where Qij ^ K''' - K'^'. If q,-/ has a mag- 
nitude of the order of the radius of the Brillouin zone, as is likely in 
most cases, the energy ftcoij' of this phonon will be a major fraction of 
of A;0, where is the Debye temperature. This is usually ^fcT in the 
extrinsic range. One must therefore use the Planck, rather than the 
Ray leigh- Jeans, distribution function for these phonons. At very low 
temperatures, inter-valley scattering is negligible: absorption of an ij 
phonon is rare because few such phonons are present; emission is com- 

' See, for example, W. Shockley, Electrons and Holes in Semiconductors, (Van 
Nostrand 1951) p. 276. 
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parably rare because few carriers have energy enough to create such 
a phonou. With rising temperature inter-valley scattering becomes 
more important. This causes t (hence n) to decrease more rapidly with 
increasing T than it would if there were no inter-valley scattering. In 
this section we shall develop this idea quantitatively. 

The matrix element for scattering of a carrier from some state in 
valley i to another state in valley;, by absorption or emission of a phonon 
hic, has the form common to all one-phonon scattering processes^ 

^■■>=fr+,v4^«f°^ ("''"'"" (18) 

(A' + 1) J (^emission *■ •' 

where A'' is the number of phonons of the given type present in the 
initial state and where Da is independent of the occupation of the phonon 
states. For inter-valley scattering D.-y is practically independent of the 
locations of the initial and final states in their respective valleys. In 
general, of course, /),> will be different for the different branches of 
the vibrational spectrum. The transition probability from a state of 
energy e in valley i to a state in valley ,/ of energy e' ^ t -j- hoi (absorp- 
tion) or e - km (emission) is proportional to | Tl/.y f times the density 
of states at energy e' in the jth valley, provided the variation of hu 
with position in the valley is negligible, as is the case for most transitions. 
Since the number of states between e' and e' + de' is proportional to 
Ae'^'^ dAe', where Ac' is the distance from the band edge, this transition 
probability has the form 

absorption: Wa == n — rr^ fl7l 

exp (hio/kT) - 1 ^^'^ 

emission: Wc « i for Ae > hin 

1 — exp (— fiLo/kT) ng\ 

for Ae < hu 

Since either of the processes (17) and (18) randomizes the initial 
velocity of the charge carriers, and since in this paper we are assuming 
the existence of an effective relaxation time r.v(e) for randomization of 
velocity by the intra-valley scattering of acoustic modes, the total 
relaxation time for lattice scattering is given by 



1 1 ' 

- = — + E [TF„(y, «) + W.(ij,a)] 



(19) 



where a labels the branches of the vibrational spectrum and W^ and W» 
' See, for example, Reference 7, p. 520. 
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Fig. 2 — Colli ribut.ions to the reciprocal relaxation^time of a charge carrier, due 




There are many transitions from valley 1 to dilTcrent ones of the other valleys, 
due to different branches a of the phonon spectrum. 

are given for each type of transition by (17) and (18) respectively, with 
ftct) = hoi(ij, a). The prime on the summation means that when a is an 
acoustic branch, the term ./ = i is to lie omitted, tlowever, since (17) 
and (18) apply to intra-valley scattering Ijy modes of the optical branches, 
such scattering is included in (19) as the terms with j = i. Fig. 2(a) 
shows the various contributions to 1/t as functions of the initial energy 
Ac of the carrier Ijeing scattered: I/t,-,- is proportional to At , as in the 
simple theory (this corresponds to a mean free path independent of 
energy for any given direction of motion), and each of the other terms 
is proportional to some (Ae ± ftoi) ". 

We shall try to estimate the order of magnitude of the steepness of 
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the parabolas describing the various inter-valley tenns, relative to that 
of the intra-valley tei-m I/tj, . For the low-frequency acoustic modes 
involved in intra-valley scattering the factor D,-y in the matilx element 
(10) is proportional to q/(ho)f''\ and since to '^ q and TV w kT/hcc » 1 
for such modes | M-a f ^ T and is independent of q. The g's involved in 
inter-valley scattering will usually be too large to satisfy kT/ko} » 1 , 
at least in the extrinsic ranges of Ge and Si, but we may hope to estimate 
a plausible order of magnitude for their WaS and WJs by assuming 
their Di/fi to be ^q/hu) with a factor of proportionality of the same 
order as for intra-valley scattering. With this assumption the steepness 
of a typical {ij) parabola corresponding to phonon emission {W^) should 
be of the same ordei- as the steepness of the intra-valley parabola when 
kT ^ hoi{ij), while for kT < hw{ij) the We{ij) parabola should become 
nearly independent of T, as contrasted with l/T,-,-{e) cc T. The parabolas 
corresponding to phonon absorption are of course always less steep, the 
ratio of the steepness of W^{ij) to that of W<.{ij) being 
WJ{U + h'xY^'' 1 - exp (- h<^/kT) 



TF./(Ae - h<^) 



i/a — 



exp {}m/kT) - 1 



= oxp{- hoi/kT) (20) 
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fie as Fig. 2, but for the simplified model of Equation (24), on 
rical calculations of Figs. 4, 5, 8, 9 and 10 are based. 
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Because of the large number of terms WaAv, «) '^ (19) °^ Fig. 2, 
each with an at present unknown amplitude and critical frequency, it 
would be pointless to undertake calculations taking individual account 
of aU the possible types of transitions. However, it is reasonable to hope 
that the behavior of the inter-valley terms can be roughly approximated 
by a model which considers absorption and emission of just a single type 
of inter-valley phonon. This model is illustrated in Fig. 3. It contains 
three adjustable parameters Wi , w^ , and hu, defined by 



/Ae\" /kT\ 



w. = 



We = 



exp (hc^/kT) - 1 



^^ (^ " ^ 



1/2 



1 - exp (- hWkT) 



orO 



(21) 



(22) 



(23) 



Equation (19) becomes 



WlT = 



At 

tlCil 



1/2 



AtV^'/kT\ -m 

ho:) \kw ) wi Lexp {hui/kT) - 1 



Ae 

hw 



1/2 



orO 



1 - exp (- ho}/kT)J 



(24) 



Thus WiT is a function of the two variables Ae/^w and kT/ho), and the 
single parameter w^/wi . The behavior of the mobility as a function if 
kT/ho} therefore depends, apart from the constant scale factor wi , only 
on W2/W1 . 

Fig. 4 shows the results of some calculations of this mobility-tempera- 
ture relation, made by numerical evaluation of (24) and (14). It is evident 
that with reasonable values of Wh/wi , the negative exponent describing 
the temperature variation of the mobility can be increased to a value 
considerably above the % of the simple theory, over a considerable 
range of temperature. This is often what is needed to explain the ob- 
served mobility behavior. In Sections 4, 6, and 8 we shall see the extent 
to which this mobility exponent is correlated with, respectively, the 
electronic part of the thermoelectric power, the ratio of Hall to drift 
mobility, and the magnitude of the magnetoresistance. 
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6 8 10 



FIr. 4 — Mobility-temperature curves for pure lattice scattering, aa obtained 
from the simplified expression (24) for the relaxation time. The quantities «;, 
and j(M measure the strength of the coupling of the carriers to intra- and inter- 
valley modes, respectively; w is the frctiueney of the inter-valley mode. The curves 
have been drawn so as to smooth out irregularities in the computed points. 
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4. THERMOELECTRIC POWER 

The thermoelectric power of a semiconductor is a little different from 
the other effects discussed in this paper, in that it involves not only the 
response of the distribution function of the charge carriers to a per- 
turbing temperature gradient or electric field, but also the alteration 
of the distribution fimction of the phonon system.''' '" The duality mani- 
fests itself in the apperaance of two contributions to the thermoelectric 
power Q: the measured Q is the sum of an electronic part Q^ , representing 
the emf necessary to counteract the tendency of charge carriers to diffuse 
from hot regions to (^old, and a phonon part Q,. , representing the emf nec- 
essary to counteract the drag exerted on the carriers by the phonons 
which flow down the temperature gradient in thermal conduction. As 
the present paper is devoted to effects having to do with the response of 
the electronic distribution function to various influences, and as all 
aspects of the theory of thermoelectric power have been discussed else- 
where,'" we shall limit the present section to a discussion only of the 
electronic part Q^ , which, fortunately, predominates greatly over Qp 
at high temperatures. 

The expression for Qr is most simply derived by making use of the Kel- 
vin relation Q, = U,/T between Qe and the eiectnuiic contribution 11, to 
the Peltier coefficient, which represents the energy flux, relative to the 
Fermi level, whi(^h accompanies the transport of unit charge in an 
isothermal conduction process. For an intrinsic semiconductor with low 
carrier concentration 

rTQ, = cn« = €r - £(, - Aer (25) 

where as before €p is the Fermi le\'el, €6 the band edge energy, and where 
AtT is the average energy of the transported electrons relative to the 
band edge, a tiuantity >0 for n-type material, <0 for p-type, and of the 
order of magnitude of kT. Now | cf - tb | can be expressed in terms of 
the carrier concentration n and the effective masses. For a many-valley 
model the luimher of (carriers n''' in each valley is easily shown to be the 
same as for a simple model semiconductor with the same | €f — ei. 1 and 
with an effective mass equal to the geometric mean of the principal 
masses mi*, ma*, win* of the valley. This is because the density of states 
in energy is proportional to the volume of K-space inside an energy sur- 
face, a quantity which for a spherical surface goes as the cube of the 
radius, and for an ellipsoidal one as the pi'oduct of the principal semi- 

» H. P. R. FreckM-ikse. Plivs. Rfv., 92, p. 248, 1953. 
"• C. Herring, Phys. Rev., 96, p. 1163, 1954. 
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axes. The total carrier concentration is therefore 

\\here Nv is the number of valleys. The final expression for Q« obtained 
bj' expressing 1 t>- — tj | in terms of n and inserting in (25) is 



8G.2 \fn m>i^ + 3^ ^n Ny + H ^n f"^.^.^") 
L H \ m m m / 



+ ^^nr+l|^lM^/deg. 



(26) 



where n is in em ^ and the upper sign is for n-type material, the lower 
for p-type. 

If a relaxation time exists, dependent only on energy, the distribution 
function for isothermal conduction has the form/**" + E-f'^' worked out 
in Section 2, and we have, for a cubic crystal, 

_ jAeVf'^'rfP ^ <^e-r> 



^'r I - .__,a) ,^ - - TTTTT (27) 



by (5), ((i) and (10), where as before the angular brackets deuote Maxwel- 
lian averages as defined by (13). 

It is important to know the value of (27) as accurately as possible, 
in the temperature range where Qe is measurable, since if (27) is known 
the measured Qc can l)e used with (26) to give information on the effective 
masses. For pure intra-valley scattering, (27) has the value 2kT. Im- 
purity scattering increases Aer by causing the current to be carried more 
by fast carriers and le.ss by slow; inter-valley scattering has the reverse 
effect. It is worth while to try to correlate the effect of inter-valley 
scattering on Atr with its effects on two measurable properties, namely, 
the temperature variation of mobility (Section 3) and the ratio of Hall 
to drift mobility (Section 6). Accordingly, calculations of | Aey | have 
been made using the expression (24) (model of Fig. 3) for the relaxation 
time. The results are shown in Fig. 5, which shows | Aer \/kT as a func- 
tion of kT/hta. 

5. PIEZORESISTANCE 

As we have just seen in Section 2, the quantum states in any small 
region of wave inimber space make a contribution to the conductivity 
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Fig. 5 _ Values of the thermoelectric transport ratio | Aer \/k,T defined by (27) , 
for the simplified lattice scattering law (24) . The ratio liij/wii measures the strength 
of the coupling of the carriers to inter-valley modes in terms of that for intra- 
valley scattering; u is the frequency of the inter-valley modes. The curves have 
been drawn to smooth out irregularities (severe for kT/ho) = 1 to 2) in the cal- 
culated points. 

which depends on (i) the degree to which these states are populated in the 
equilibrium distribution/"", (ii) their group velocity, and (iii) their relaxa- 
tion time, or more generally, the transition probabilities for scattering 
from these states to others. When the crystal is strained, any or all of 
these factors may change, and the resulting change in the sum of all the 
local contributions to the conductivity constitutes the piezoresistance 
effect recently discovered by Smith.'' Although there are a number of 
processes which can contribute to the three factors (i) to (iii) just 
enumerated, it can be argued plausibly that for a simple many-valley 
model the principal effects are usually those due to a single process, 
namely, the strain-induced shifts of the energies €»■ of the band edge 
points K'*'. We shall consider (i) to (iii) in turn: 

(i) The change in the population /*''^(K''' + AK) depends on the shift 
of the energy ((iv'*' + AK), and because. of the smallness of AK this is 
practically the same as the shift 5e''' of e(K'''). In a shearing strain some 
of the 5e''' will be positive, some negative, and so some of the valleys will 
have their populations decreased, some increased, the fractional change 
in each case being be'^/kT. Now it is evident from the second equation 
of (7) that the contribution of a single valley to the conductivity is aniso- 
tropic. If all valleys are populated equally, as we assumed in Section 1, 
the total conductivity will be isotropic. But if strain causes different 
valleys to have different populations, the overall conductivity will have 
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an aiiisotropy like that of the more populous valleys. If the ratio of 5e''^ 
to the shear strain amplitude is of the order of magnitude of the known 
ratio of Sea to st,riiin for isotropic compression, viz., a few volts, the frac- 
tional change of/"" per unit strain will be of the order of hundreds. Since 
the observed fractional change in resistance per unit shear strain is of 
the order of 10 in the more favorable orientations,^ the change in/"" is of 
the right order of magnitude to contribute a major part of the effect. 

For an isotropic compression or dilation there exists the possibility, 
not present for shearing strains, that the total carrier concentration may 
be changed in first order. A large effect of this sort, again of the order of 
the 5t^'^/kT, will occur for a specimen in or near the intrinsic range, be- 
cause of the change of energy gap cg mth strain. This effect rapidly 
become.s negligible, however, as the speehnen is made extrinsic. For ex- 
ample, if unit volume of an n-type specimen has an excess no of donors 
over acceptors, all ionized, the hole and electron concentrations m, , 
Ue , obey 

lit, = nu + rth , rieni, = nl 

where n,{r) is the value of Ue = iik in intrinsic material. Thus if Ud » ?i, 



. , (28) 

and since dnjUui = dn-h/dea , the energy gap effect is negligible if iiu 
exceeds ni by a large factor, even though the change in «,- with strain 
may be sizable. For extrinsic specimens with incomplete ionization of 
impurity centers, there may of course be an effect of compression on total 
carrier concentration due to change in the ionization energy of the 
centers; however, if this ionization energy is <^€g this effect will be of a 
smaller order of magnitude than the he'^ fkT. 

(ii) It is easy to show that the fractional change in group velocity 
per unit strain must be much smaller than the Se'^/kT just discussed, 
hence too small to contribute in a major way to the piezoresistance ef- 
fect. For we expect the change 5v in the group velocity at K''' + AK to 
have an order of magnitude given by 

Bv — [5e(K^" + AK) - 5e^'^]/hAK — {AK/K^'y(St^'^/hAK) 

since the quantity in square brackets must vary as AZv^ Since 

V ^ [€(K"' + AK) - t^'^]/hAK = U/hAK 

we have 

Bv/v ~ (AA7/v^'')'(5e'''VA0 (29) 
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Since a typical charge carrier has Ae ^ kT, (29) is smaller than the ratio 
discussed in the preceding paragraph by the factor (A/v/A^"')'. It is thus 
plausible to neglect strain-induced changes in group velocity, or equiva- 
lently, in the effective masses. 

(iii) Consider the transition probability from a state K to the group 
of states lying in a small element of volume in K-space, centered on a 
point K' at which the proper energy conservation law for the transition 
K -^ K' is satisfied. This probability, like all quantum-mechanical 
transition prolmbilities, can be expressed as the product of the square of 
a matrix element M{K, K') by the number of states per unit energy in the 
given element of volume. We have to consider the effect of strain on each 
of these factors. 

The matrix element M(K, K') can be changed either by a change in 
the wave functions -^k , *k' , or by a change in the physical processes 
deteimining the perturbation operator M, e.g., a change in the ampli- 
tudes of the thermal vibrations, or a change in the dielectric constant, 
which enters into scattering by charged impurities. Typical assumptions 
on the e(|Uatiou of state of a crystal suggest that the fractional change 
in the squared vibration amplitude, per unit strain, might be of the 
order of a few units, i.e., at least an order of magnitude less than the 
observed elastoresistance for the optimum orientations. The effect of the 
change in the wave functions is of similar magnitude: To effect a major 
change in ilf (K, K') one must make a major change in the wave func- 
tions. To do this probably usually reciuires a strain of amplitude 0.1. to 1 . 
Therefore it is reasonable to expect that the fractional change in | Al' | 
per unit strain will be of the order of 10 or less, i.e., again an order of 
magnitude smaller than bt'^/k'l\ or than the observed elastoresistance. 

The effect of strain on the density-of-states factor, on the other hand, 
can be larger. For intra-valley scattering, where initial and final states 
are both near the same band edge point K*", the effect is of course very 
small, since initial and final states undergo very nearly the same energy 
shift with strain. But for scattering from one valley i to another valley 
j, the two energy shifts 5e''' and 5e*'' are in general quite different, and 
for a given initial state apphcation of a strain will change the set of K"s 
describing final states which conserve energy and hence will change the 
density of final states — e.g., the density in a given solid angle of vectors 
AK' = K' — K^^'. Since in a given solid angle the density of states is 
o: Ae'""', the fractioiial change in this density due to a strain is Se'^'/SAe', 
which on the average is of the order of he'^ /kT, i.e., of the same order as 
the effect discussed under (i). 
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Table I — Ways in Which 


Sth- 


\iN Can Affect Conductivity 


Effect 


Probable Order of Magnitude 


Rank in 
Importance 


(i) Population function f-"^ 




U^'>/kT + much 
smaller terms 


First 


(ii) Clroup velocities of states. . . 




«6t<'VA-r 




(iii) Tninsition probabilities 








(a) Matrix elements 








(a) Wave functions .... 




^5e<'-VA-T 




09) Viliration amplitudes, 


etc. . 


Rather <6ci'^/kT 


Second (?) 


(b) Density of states 








(a) Intravalley 




«5t">/A-r 




(/3) Inicrvalloy 




5e("A7' + much 
smaller terms 


First 



Table I summarizes the foregoing discussion of the ways in which 
strain can alfeet conductivity. 

Appendix C gives the mathematical treatment of the two effects which 
are of the order of the quantities Se^^^/kT, namely, the change in/"' 
and the change in the density-of-states factor in the transition proba- 
bilities for inter-valley scattei'ing. This treatment, which is fairly simple 
and straightforward, is based on the following assumptions: 

(a) Neglect of all other effects of strain on the conductivity. 

(b) Tho assumption of the preceding sections that the scattering of 
tlic carriers is describable by a relaxation time which in each valley is a 
function of energy only. 

(c) Carrier concentrations in the extrinsic range. 

(d) Maxwell-Boltzmann statistics. 

(e) Valleys lying along a threefold or fourfold symmetry axis of a 
cubic crystal. For such valleys the energy surfaces are ellipsoids of revo- 
lution. 

The principal features of the calculation are qualtative ones which 
can be derived with little or no mathematics. These we shall consider 
here, with a little inquiry in each case as to the sensitivity of the conclu- 
sion to relaxation of the assumptions (a) and (e) above. The first such 
feature to be noted is that under assumption (a) the. change of mobiliti/ 
in an isotropic compression vanishes. For in an isotropic compression all 
the liand edge shifts 8e '' are equal. This means that for a given total 
carrier density the distribution function /*"* in each \'alley does not 
change, if, as we are doing, we neglect changes in the effective masses. 
Similarly, since all valleys are shifted together, there is no change in the 
density of final states corresponding to any inter-valley scattering 
process. The present conclusion is easily seen to be independent of 
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Table II — Isothermal Elastoresistance Constants for Ge and 
Si (Smith, Reference 4) 









I da Hii[ + 2hiii 


Material and Resistivity 


fM12L! ^ »l1t 


2 ~ 2 


ndlnV 3 


licm 








nGe 1.5 


-93.0 


+0.4 


-5.3 


5.7 


-92.0 


+0.5 


-6.8 


9.9 


-92.8 


+0.1 


-9.8 


16.6 


-93,4 


+0.1 


-13.6 


p Ge 1.1 


+65.1 


-2.8 


+3.9 


15.0 


+66.5 


-6,3 


+ 1.4 


n Si 11.7 


-10,8 


-79,5 


+5.7 


p Si 7.8 


+ 110.0 


+3.9 


+6.0 



Here m„,^p , defined by (C7) of Appendix C, describes the relative change of the 
conductivity tensor with the strain tensor, in a coordinate system oriented along 
the cube axes. The abbreviation of this by m„{r, s = 1 to 6) follows the same prac- 
tice as that used for elastic constants. 



assumptions (b), (d), and (e). As regards assuniption (a), however, it is 
clear that inclusion of any of the other strain effects listed in Table I will 
in general lead to a nonvanishing effect of compression on the mobility. 

By virtue of the fact just mentioned it is possible to test the vahdity 
of assumption (a) by comparing the observed elastoresistance for iso- 
tropic compression with that for a typical shear. Table II, taken from 
the work of Smith,'' shows the room temperature elastoresistance con- 
stants of Ge and Si. The entries in the last column vary with resistivity 
for the case of Ge, because of the energy gap effect discussed under (i) 
above [failure of assumption (c)]; our present interest is therefore in 
the values for low resistivity specimens. For these the volume coefficient 
(last column) is in all cases only a few percent of the larger of the shear 
coefficients (middle columns); this accords with the expectation that 
the volume variation of the squared matrix element for scattering (pre- 
sumably the largest of the neglected effects) should be an order of magni- 
tude or more smaller than the de^'''/kT. This is encouraging, but it must 
be remembered that the shear variation of the matrix element may well 
be larger than its volume variation because suitable shearing strains can 
usually couple a band edge state to states closer to it in energy than can 
isotropic dilatation. 

The second important conclusion is that tinder assumption (a) the 
change of mobilitij vanishes for a dilatation along a (1 00) direction if the 
valleys are on (/?7) axes, and for a dilatation along a {111) direction, if the 
valln/s are on {100) axes. In terms of the elastoresistance coefficients of 
Table II, (mn — mi^) vanishes for (111) valleys, and tuaa vanishes for 
(100) valleys. This conclusion is obvious from the symmetry of the 
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problem: a shear compounded out of a unidirectional dilatation of the 
type described and an isotropic compression must shift all band edge 
points by the same amount. This amount must be the same as in the 
negative of this shear, so all 5e,- = 0. This conclusion is again independent 
of assumptions (b) and (d), but in general breaks down if assumption (a) 
is relaxed to the extent of taking account of the effect of strain on the 
matrix e'ement for scattering. 

The third point to be made is that the change of mobility accompanying 
a given strain is inversely proportional to T at temperatures low enough for 
inter-valley scattering to he of negligible importance. This is because the 
relative change of population of different valleys with strain is propor- 
tional to the Se-'^/kT. The more complete treatment of Appendix C shows 
that, under the present assumptions, the decrease of elastoresistance with 
increasing T should be more rapid than 1/7* when inter-valley scattering 
is just becoming important, but that for very high T it should again go 
as l/T. This behavior is illustrated schematically in Fig. 6. The present 
conclusion is not dependent on assumptions (b) or (e), but depends on the 
others, especially (a). The effect of strain on the matrix elements for 
scattering will give a contribution to the elastoresistance which is inde- 
pendent of T in the range (if such exists) where only intra-valley lattice 
scattering is important; if impurity or inter-valley scattering contributes 
the dependence is of course more complicated. 




LOG T 



Fig. 6 — Schematic variation of any component of elastoresistance with tem- 
perature, showing the transition from the low temperature region where the only 
important effect of strain is to change the relative population of the valleys, to 
the liigh temperature region where the effect of strain on inter-vailey scattering 
is of comparable importance. 



258 THE BFXL SYSTEM TECHNICAL JOURNAL, MARCH 19o5 

In conclusion, a few words are in order regarding the extent to which 
the present conclusions on piezoresistanee can be expected to hold for 
models other than the simple many-^-alley type to wliich this memoran- 
dum is restricted. First of all, we may note that for the various types of 
degenerate models, such as that of Fig. 1(a), different states inthe neigh- 
borhood of the same band edge point can experience widely different 
energy shifts 5e under shear. Consequently the mobility will be affected 
in a major way not only by the two effects labeled "first" in Table I — 
change of the population fmiction and change of the density-of-states 
function in inter-valley scattering — but also by changes in the group 
velocities. Moreover, the perturbation of the crystal Hamiltonian by a 
uniform shearing strain may now have sizable matiix elements between 
states of the same wave vector belonging to the different bauds which 
come together at the band edge point. This can cause the form of the 
crystal wave functions to be much more sensitive to strain than when the 
perturbation only connects states a few volts apart in energy, and so 
the dependence on strain of the matrix element for scattering may be 
much larger than for the simple many-valley case. Thus at least four, 
rather than two, of the entries in Table 1 become of first magnitude. 

In view of these facts, most of the conclusions reached for the simple 
many-valley case probably become invalid for the degenerate and de- 
generate many-valley cases. An exception is the conclusion concerning 
the smallness of the change of mobility in isotropic compression. The 
perturbation introduced into the crystal Hamiltonian by an isotropic 
compression does not mix states of a degenerate set, so the arguments 
previously given remain valid. 

6. HALL EFFECT AT LOW H 

When a magnetic field is present the E in the transport equation (3) 
must l)e replaced by E + v X H/r, where as before, v is the group veloc- 
ity. Thus with the upper sign for electrons, the lower for holes, the dis- 
tribution fmiction / of the carriers obeys 



■' ^ (30) 



dt [_ c J T 

We shall seek the solution of this as far as the first order in E and the first 
order in H, i.e., we shall set 

/ = /<«' + E.f'^"'+E^V^J./'"+ . . ■ (31) 

IIP 

There is, of course, no term of the first order in H and the zeroth in E, 
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since a pure magnetic field has no effect on/. The vector function f '" is 
of course the same as in Section 2, Equation (5), namely, 

f^'"' = ^erV^f (32) 

Putting (32) and (31) into (30) we get for/^/"^ 



c 



whence 



r (11) _ cV -^ d {df^°^\ 



(33) 



where S^„0 = if any tuo of its suffixes arc the same, and ±1 if the suf- 
fixes are an even (odd) permutation of xyz. 

The physical meaning of the steps leading to (33) is just that a weak 
magnetic field perturbs the /'" solution of Section 2 by displacing each 
part of the distribution in the direction of v X H in crystal momentum 
space, the displacement being proportional to v X H and to r. 

The term (33) in the distribution function gives rise to a contribution 
j*'" to the current, which is at right angles to H and to E. This contribu- 
tion can be described by a "Hall conductivity tensor" ff>^^ , thus: 

J\''" - E ^x..B.//. (34) 

The contribution ax^/" of the Vth valley to o-.^, is easily obtained from 
(33). We shall assume Maxwellian statistics, so that df^^/dP^ = 
-{vJkT)j"^\ When this is inserted into (33) the last factor involves a 
derivative of vj^^\ with respect to P„ . If r depends only on energy, as 
we are assuming throughout this memorandum, the derivative of / t 
with respect to l\ in (33) will be proportional to v„ , and Sa.s^.as''^''" 
will vanish identically because of the anti-symmetry of 5,„fl in a and IS. 
Therefore the only term which need be retained in d^dP^ is that in 
dvJdPa . If the coordinate axe.s are chosen along the principal axes of 
the energy surfaces of the iX\\ valley, this latter derivative is just h^Jm^*. 
Thus we get, with the upper sign for n-type the lower for p, 

3 1 / /(O' 



OhJ'^ = =F- E ^ X! K<^aVKV0 



f-"'.. 
V "» 



Z^ ' i" ^ ""fi^ Tt^ 



(35) 
where as usual the first summation is over all vectors AP^" in the )'th 
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valley, per unit volume, and over both states of spin. In our present co- 
ordinate system the average of V\Vff over an energy shell vanishes unless 
^ = X^ while that of vk can be evaluated from the equipartition relation 
(10); v\ — » 2Ae/3mx*, where Ae = | e — es [ is the distance from the band 
edge. Thus with kT = H <Ae>, (35) reduces to 

,,,."> = T i!^ %L> 4-_ (36) 

c (Ac) mx*in^* 

where n^*', as in (9), is the number of carriers in the ith valley per unit 
volume, and where the angular brackets are Maxwellian averages, as in 
(13). 

The propoi'tionality of the Hall condxictivity tensor to <AfT > and to 
the reciprocal product of two different principal masses is easy to under- 
stand physically. Without a magnetic field, an electric field in the fx 
direction gives a distribution, in each energy shell of the ith valley, 
which has a mean velocity in the fi direction proportional to Aer/m^* 
(of. Section 2). Thus the distribution in this energy shell is acted on by a 
transverse magnetic force whose average value is proportional to this 
expression. This transverse magnetic force produces a transverse current 
proportional to the force and to r/mx*, where X is the transverse direction. 

For a cubic crystal the relation of the Hal! current to E and H must 
Ijc isotropic, i.e., the right of (34) must be proportional to E X H. It is 
easily shown that the quantities in (34) are related to the ordinary con- 
ductivity (To , Hall coefficient R, and Hall mobility hh = TtltToc, by 

j^"' = a,'R E X H 

or 



o'^M•' 



= tro RSx^v — -F B\pr (37) 



where as usual the upper sign is for n-type, the lower for p. Since 
S),^,.o-j,^^5j,,.^ is invariant with respect to changes in the orientation of the 
coordinate system, we may evaluate it by evaluating each S>^^ff|^f ' 5),^^ 
in the system of principal axes of the ith valley, and then summing on i. 
From (36) we find in this way 

tro 72 = =F = Vq Z^ a■^M^5^;-F = ^--i A^ 2-, fAc" ^Ajiw 

/" , ' '" , (38) 

en <AeT > ^ / 1 , 1 , 1 



c <Ae> ^ \mi*m2* m2*in3* m*7ni 
where n = S.^i''' is the total density of carriers. A neater way of pre- 
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Table III — Values of the Last Factor in (39), for Cases of the 

Form 



vii* ^ ?»2* = »u*, in* = V),,*. 









20 

10 

5 

3 

2 



0.784 
0.816 
0.868 
0.918 

0.960 



1 


1.000 


0.5 


0.938 


0.3 


0.808 


0.2 


0.674 


0.1 


0.437 


0.05 


0.2M 



senting this result is in terms of the ratio mh/m- Multiplying (38) by c/ctom 
and using ao = ne^ and (14) for /i we get 



/iff 


= 


<AeT'><Ae> 


fi 


<Aer>- 






<A€rS<Ae> 



1 



niim-i 



+ 



1 



ni2 m-A 



+ 



[m*'^^'^^) (39) 



<Acr>= 



■B, 



say. 



Note that the first factor of (39) is the value of mh/ju in the simple 
theory/^ and that the second factor B, involving the ani.sotropy of the 
elfective mass, is unity for zero anisfjtropy and < 1 in general. Some sam- 
ple ^'akies of this mass factor B are given in Tabic III and Fig. 7. Fig. 8 
gives values of the first factor in (39), for the simplified model of Intra- 
and inter-valley scattering described by (24). 

7. THE BENEDICT-SHOCKLEY EXPERIMENT 

We turn now to the response of tiie assembly of carriers to an electric 
field which varies sinusoidally with time. As Benedict and Shockley have 
shown, this response becomes limited at high frequencies by the inertia 
of the carriers, and so by measuring it one can obtain an cffecti\'e mass. 

" See, for example, Reference 7, |>. 277. 

" T. S. Benedict and W. Shockley, Phys. Rev., 89, p. 1152, 1953. 
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Fig. 7 — Dependence of 



B = 



mi"wi2 



+ 



mi'iUs 



+ 



m^'mi 



I \ 1 1 \^ 

( ^ + — * + — J 

on the anisotropy of the effective mass, for the case »ii* = "'2* = "U*> "'a* = »'*■ 

The solution of the transport equation for this case proceeds almost 
exactly as in Section 2. We assume that the scattering of the carrier is 
described by a relaxation time t, whose dependence on position in mo- 
mentum space we shall for the moment leave unrestricted. The analysis 
starts as before from (3) for the distribution function / of the carriers, 
namely, with the upper sign for electi'ons, the lower for holes. 



at 



(40) 



(We neglect the very small effect of the magnetic field generated by 
dE/dt.) Instead of (4) we write, if E = Eoc'"', 

/(O = /"» + i'o-f'"(0 + O(En') (41) 

From (40) and (41) the equation for/'^' is 



dt 



r 



T 



(42) 
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Ah 




Fig. 8 — Values of the ratio 






<Af7->2 



for tlie aimplified lattice scattering law (24). The ratio w-Jwi measures the strength 
of the coupling of the carriers to intcr-vailey modos in terms of that for intra- 
valley scattering; w is the freijuency of the inter-valley modes. The curves have 
been drawn to smooth out irregularities {severe for kT/hu = 1 to 2) in the calcu- 
lated points. 

If Ae = I t ~ tft I is the distance from the band edge, we have, for Max- 

wellian statistics, /*"' cc exp (—Ae/A'T) and of course V^Ae = v,thegroup 
velocity. Thus f*" = /o*"e'"' with 



In 



kT{i + ZW7-) 



(43) 



The current density is given by the usual sum. over the different valleys 
i and the states (AP''\ spin) in each valley: 



i = E E (=Fe)v(AP"')/(AP^'') 



From (41), (43) and (44), 



(44) 



(45) 



7^ jp('"),,, \/;7'(l + i^T 

3/;7'r ap4t',,„V(1 +^W)/' 
if the crystal has cubic symmetry. Thus the semiconductor has the fre- 
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quency-dependent complex conductivity 



«« / ^-._\ (46) 



where n is the number of carriers per unit volume and the angular 
brackets denote a Maxw^ellian average. Note that in deriving (46) we 
have not had to assume anything about the dependence of e or r on 
AP; (46) is therefore valid for all models, not merely for the many- valley 
case. However, (46) is still not exphcit enough to be directly usable for 
the evaluation of experimental results, and we shall need to use the spe- 
cial properties of the many -valley model to express the Maxwellian aver- 
age in terms of measurable or readily interpretable quantities. 

Under the usual assumptions of this paper r is a function of energy e 
only, so the average in (46) depends only on the average of v over an 
energy shell. By the equipartition principle (10) we may therefore replace 
u^by 

^ / W»/ + -A^ + ^A^\ = y^^, (^ + _L + J,) (47) 
\ mi* ma* ma* / \mi* ma* mg*/ 

The average of the masses is the same m'^' we encountered in Section 2, 
namely, 

— = Mi—.'r—^^^ (48) 

m"^' \mi* m-a* ma*/ 

Therefore 

/ .V \_ 2 / A.r \ 
\l + lojr/ m'^' \l + zcor/ 

The expression for a{ui) becomes, with its real and imaginary parts sep- 
arated. 



, , Ine \ I Aer \ . I Act \ 



(50) 



The real part trK(w) of this is what is usually called the "conductivity". 
The imaginary part ai{ui) is proportional to a contribution to the dielec- 
tric constant /f(aj), since (4T)~'K(a))6E/a/ is the sum of the displacement 
current (47r)~VoaE/3i and the part of the true current j' which is in phase 
with SE/dt. Thus the departure of h-(w) from the dielectric constant ko 
of the crystal without its free carriers is given l^y 



KO 



— K{(ii) = —- — 0"/(w) (51) 
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At the frequencies and temperatures which have been used for the 
Benedict-Shockley experiment, most of the carriers have relaxation times 
short compared to iiT , and it is appropriate to make an expansion in 
powers of w: 

( l +'LvV " ^^'''^ " c,\Aer'> ■ ■ • (52) 

(If T ^ CO as AP — > 0, as would be the case if the only scattering were by 
phonons of negligible energy, the series (52), (53) do not converge for any 
finite w. However, asymptotic series can be written down which differ 
only in order u and higher from the series obtained by simply expanding 
the denominators) . Denoting the dc conductivity by co — equal to ne 
times the ju of (14) — we have from (50) to (53) 



.«W=..[l-^^+0(c.')] 

Ko — li{lJl) = 4x0-0 ~— — 0} — + O{0) ) 

\_ <AtT> <AeT> 



(54) 



(55) 



It is convenient to express the first tci-in in the sqiiare bracket in (55) 
in terms of the Hall mobility /xh , since the same average <AeT^> occurs in 
(38) as in (55). Let us set o-q = nc^i and use the designation B for the 
last factor in (39), a factor ^ 1 dependent on the anisotropy of the 
effective mass in each valley and close to unity unless the anisotropy is 
very extreme (see Table III). Then 



Kg — k{u) = ^ 



1-- -^B- <Aer^>^ +^^"^J f^^^ 



Equations (55) and (56), like all equations in this memorandum, is in 
Gaussian units. For rationaUzed MKS imits, as used in the papers of 
Benedict and Shockley, the coefficient 4it should be replaced by 1/eo , 
where eo is the permittivity of the vacuum. 

The leading tei-m of (56) is the same as that which one would obtain 
by simply replacing /i" by {ihh/B in the formula used by Benedict and 
Shockley (simple model, t ^ constant). But because the dimensionless 
factor <AeTS<A€>/<A€r'>" is always ^ 1 instead of ^ 1, the second term 
in the bi'ackets in (56) is not the same as that resulting from this substitu- 
tion. Thus the expre.ssion used by Benedict" in his later analysis of data 

'3T. S. Benedict, Phys. Rev., 91, p. 1565, 1953. 
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Fig. 9 



Sample values of the quantity ^^^ ^^, occurring in the equation 



{56} for the high-frequency dielectric constant, for the simplified hitticc scatter- 
ing law (24). The ratio xv-Jivi measures the strength of the coupling nf the carriers 
to inter-valley modes in terms of that for intra-vidley scattering; w is the fre- 
quency of the inter-vallcy modes. The dotted extrapolations are qualitative only. 

Oil p germanium is corrected only to the zeroth order in oj and to the ap- 
proximation .^ f^ 1. Some sample values of the ratio <AeT ><Ae>/<AeT > 
have lieen computed for the scattering law (24), and are graphed in Fig. 
9. These show that the range of possible variation of this factor is con- 
siderable. 

A similar, though less useful, transformation can be made on (54), to 
express the second term in brackets in terms of the dimcnsionh^ss cooffi- 
cieiit G defined by 



G - 



<AeT ><Ae>' 
<Aer>3 



(57) 



This is a coefficient which we shall encounter in the next section, in the 
theory of magnetoresistance, and which is graphed in Fig. 10 for the inter- 
valley scattering law (24). We find 

1 - J'^t-^G^- 0(co')J (58) 



The quantity G is ^1, the equahty holding only if r is a constant. Table 
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IV and Fig. 10 give some typical values, for scattering laws of the 
form T cc ^e or of the form (24). 

8. LOW-FIELD MAGNETORESISTAXCE 

III Section we set up the Boltzmann equation for the steady motion 
of charge carriers under the combined influence of an electric field E and 
a magnetic field H, and solved it to the first order in H. We shall now 
undertake to solve this equation to the second and higher orders in H. 
The solution has been worked out independently for a number of cases 
by Abeles and Meibnom,"' Shibuya,^ and Shockley (unpublished). We 
shall not give all the details of the solution, especially at large H, as many 
of them can lie found in the reference just mentioned. However, to em- 
phasize some features not brought out in this previously published work 
we shall review the whole calculation briefly from the beginning. 

The relation of theory and experiment in the area of magnetoresistance 
resembles that for piezoresistance, in that the tensor quantity which is 
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G = 



<A«->= 



and -I = 



<At7^><Aer> 

<A*T*>* 



for the sirnpliHcd lallk-c srat (criiiK law (24). Tlie ratio Wi/wi measures the atrerigtli 
of tlie ('oii|ilin(; nf llie carriers to iiiler-vjilley modes in terms of that for intra- 
valh'\' .seal tciiiin; tu i« llie frequciu-v of the inter-valley modes. The G curve has 
been dra\vn lo he rouKhly eonsistent with the smoother A curve and the curve of 
Fip. S. The dotted extrapolations of the curves are inteude{i only to show the ex- 
pected (|ualitative behavior. 
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Table IV — S-uhple Values of the Quantities G and A Defined 
BY (57) and (68), Respectively 



Scattering Law 




<AeT'> <i,r> 




2.58 
1.77 
1.00 
1.33 
2.52 
5.89 
2.77 
2.70 
1.92 
1.52 


4.01 


T cc At-"-^ 


1.27 




1.00 


T DC Ae"-^ 


1.09 




1.28 


T « Ae'-^ 


1.58 


Form (24), w^/wi = 4, kT/hu 


= 1.43 

0.667 

0.333 

0.167 


1.40 
1.21 
1.08 
1.15 



simplest to calculate is reciprocal to the one which is directly measured. 
Thus one measures piezoresistance but calculates elastoresi stance. Simi- 
larly the measured magnetoresistancc is a change in the electric field E 
for given current, whereas the simplest quantity to calculate is the change 
of the current for given E, i.e., the dependence of the conductivity tensor 
on H. We shall see below that the neatest way of comparing theory and 
experiment, at least for small H, is to invert the observed magnetorc- 
sistivity tensor to get the magnetoconductivity tensor, and then compare 
the latter with theory. 

The Boltzmann equation (30), from which we shall start, is 



0^^eE...f^e^-^...f-^I^I^^ 



(-59) 



where as before the upper sign is for electrons, the lower for holes, / is 
the distribution function of the carriers, /"" the (Maxwellian) distribu- 
tion in thermal equilibrium, v = Vpe the group velocity. As in Section 2 
we set 



(1) 



/=,r +E-f'> + 0(-BO 



(60) 



ed) 



and neglect higher order terms in E. The resulting equation for f can 
be written in the condensed form ' 



r(0) 



(1) 



- ±TeVi.r' ± TH-fi" -f 



eC) 



(61) 



where 



Y = - V X Vp 
c 



(62) 
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In the notation of Davis, and Seitz, y = (^/^ '^)^) while in the nota- 
tion of Abeles and Meiboom,^ Y = (e/c)U. The solution of (61) can be 
expressed formally in teiins of the reciprocal of the operator (I ± tH ■ y) - 

I f*" = ±(1 ± rH-Yr'TeVi./"' (63) 

If we set 

(1 ± tH-y)~' = 1 + tH-y+ (tH-y)(tH.y)--- (64) 

we see that the leading term of (63) is just the solution (5) of Section 2 
for // = 0, while the next is just the term (33) of Section 6. In other 
words, the series (64) corresponds to an iterative solution of (59). If we 
are interested only in the first few powers of H, this iterative solution is 
as simple as any ; at high fields it is better to solve (61) explicitly in closed 
form, a procedure we shall outline in the next section. 

The solution given by (63) and (64) of course applies for any depend- 
ence of the relaxation time t and the energy e on position P in crystal 
momentum space. However, we are here interested only in the case 
where, in each valley i, e is a quadratic function of the components of 
AP = P — P''' and where t is a function of e only. For this case some 
simplifications are possible. For one thing, t commutes with the operator 
Y, since (62) acting on a function of e contains the factor v X Vpe = 0. 
The exprcssioit for the current density j in powers of H has the form 



h 



= E <^.^F., + E <r>^>'cH.E, + E T^..^HM^E, + ■ ■ - (65) 



where of course a^y = a^hfi^ for a cubic substance, with o-o given by the 
etiuations of Section 2, and similarly c^.a = oa'RS^vc , where R is the low- 
field Hall coTistant and as in Section 6 5;,va = ±1 if nva is an even (odd) 
permutation of 123, zero otherwise. To get the contribution of the I'th 
valley to the second-order "magnetoconductivity" tensor a^^a^ , we 
multiply (60) by ±ra^ , insert (63) and (64), and sum on all momentum 
vectors in the ith valley and in unit volume, and on spins. The result is 



f^^y^e'" ^ T7i^ J^ /'''V'(Vp7a'y^y.)Bycnm (66) 



(i) ^ 

k-T i^>,. 



where as u.sual we have assumed / " to be Maxwellian, and where the 
subscript "symm" means that the expression in parentheses is to be 
averaged with the expressions obtained from it by permuting a with 0, 
since only the part of cF^y„ff symmetrical in a and ,3 has physical signifi- 

'* I,. Davis, Phys. Rev., 56, p. 93, 1939. 
■» F. Scitz, Phys. Rev., 79, p. 372, 1950. 
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cance. (This symmetrization is not necessary, but simplifies the work by 
preventing the appearance of meaningless components.) 

The cxphcit e\-ahiation of (66) is a straightforward but tedious exercise 
in algebra, and will not be given in detail here. However, there are some 
important properties of the tensor (66) which can be estabUshed rather 
simply. Since the components of v are linear functions of the components 
of AP, the operator r defined by (62) takes any linear function of the 
AP\ into another Hnear function. Therefore the v^ycy^ih in (66) is a (]uad- 
ratic function of the AP^ , and it is easily seen that this function contains 
denominators of the fourth degree in the effective masses. Xow a (|uad- 
ratic function of the AP\ can be written as an effective mas.^ times the 
energy Ae relative to the band edge, times a function of the direction of 
AP, dependent only on the ratios of the effective masses in the principal 
directions. Thus we may write, for example, 

Ae 

v^yay^f,. = "oTa X function (fiuaQ, direction of AP, mass ratios) (67) 

where m'" is the inertial average of the effective masses, defined by (15). 
Now let the summation on AP be broken up into a summation over 
values in an energy shell Ae to Ae + dAe, and a summation over different 
shells. The function of direction in (67) will be the same for all the shells, 
and so we ha"\'C the result 

<'■] — — ^ p ('■> 

where f^.^s''' depends on the anistropy ratios of the effective masses in 
the I'th valley, but not on the variation of r with energy, while a is pro- 
portional to the number of carriers in the ^'alley and to the Maxwellian 
average of t Ae. 

It is convenient to express the average of T''Ae in dimcnsionless form 
by using the (luantity G defined by (o7), namely, 



G = 



<Aer%<Ae>' 
<AeT>^ 



or else the cjuantity 



A = ^l''><^"> (68) 

<AeT"->- 



Here as usual the angular l)rackets denote Maxwellian averages us de- 
fined in connection with (14). We njay use (14) to eliminate m'" and the 
carrier density and if we wish we may eliminate fi in favor of mh liy (39). 
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Tlie result is 

where o-o , fi, n,i , are the conductivity, mobility, and Hall mobility, 
respectively, at // = 0, Nv is the number of valleys, and where B ^ \ 
is the function of the efTeetive mass ratios defined in (39) and Talilo III. 
Summing on valleys i gives 

cr,.„^ = .1 (ff<i ~) ^.-3 = G L,^ B%,.0 (70) 

/'V^otf = Xp Zj ^\yoff ' (71) 

iV I- i 

Note that F^„„^ , as defiiied by (70) or (71), is dimensiouless, as are G 
and .4 ; Fp„a^ or B' F^^as depends on the geometry of the valleys and the 
ratios of the pi'incipal masses of a valley. We shall see presently how the 
analysis of experimental data is facilitated b^' this decomposition of the 
magnetocoiiduclivity into the product of a scalar factor depending on 
the beha\'ior of t and a tensor factor depending on the shape of the 
energy surfaces. 

The quantity A dehned by (fi8), lilce G, is ^1, the e(]uality holding 
only if T is a constant. Some sample graphs of A and G arc shown in 
Fig. 10, for scattering laws of the form (24), and some numerical values 
for this case and for t cc Ae'^ are given in Table IV. Note that for the 
ideal case of intra-vallcy lattice scattering only, a case approximated 
in \'ery pure material at moderately low T, r ^ — i^ and A = 4/ir = 
1.27, (' = 9x/16 = 1.77. Table \' gives values of all the nonvanishing 
coefficients F^^^^''' i'elati\'e to a coordinate system oriented along tfie 
principal axes of a valley. The middle rows of Table \l give the f'\raff , 
relative to the crystal axes, for some of the simpler possible arrangements 
of \'alleys. The entries were olitained, of course, by comparing (()9) or 
(70) with the results of explicit evaluations of ((i6). For completeness, 
Table V also gives the dire(^tional factors involved in the contribution 
o-^„"' of a single valley to the conductivity tensor o-^,, in the alisence of a 
magnetic field, and to the Hall conductivity tensor a^^^ defined by (34) 
or (65). All these table entries are similar to those given by Abeles and 
Meiboom.'' However, they have given the unsymmetrized o-^„„3 etc. for 
the cases r = — ' ^' and +-'2, in terms of mean free path, absolute \-alues 
of the masses, and carrier concentration; here we have given the sym- 
metrized a^vaff etc. in terms of the directly observable o-o and ^ih , and for 
any 7(Ae). 



3 




































1 


^ 






^ 


" 


9 


• • 






•.% 


*H*_ 




1> 




i 




.-H 


H .- 




£ £ 
1 


::?? 


I 


£ £ 
1 


£ £ 
1 


:^ 


3 








ft 


1 






1 






oc 
































« 
























~ 




,._^ 


























# 






















• 


« 


c 












3 
1 








o 


i 
* 

+ 
« 


* 
+ 


* 

+ 
* 

» 


• 

+ 
♦ 

£ 
+ 


-1 

* 

a 

+ 


H 

e 

+ 


# _ 
+ 


g 
+ 


s 

* 


»_ 

» 

IN 

+ 


* 

s 

+ 


^ 
+ 












+ 

* 

1 


+ 


+ 




J, 


3 


CO 


s 




* 


















J 














































CO 




^ 






























^ 














"^ 








"?. 














\ 








fc,' 




























C 




a 


s 




ri 


ffi 


II 


a 








3 5 


\ 





Q 


tC' 


tC 




d 




1 


b r 


1 1 

a 
> 


to. 
to. 

a 


11 


iC 


II 


II 


II 
















• 








>^ 














H 














S 








ffH 








1 






II *^ 


a 






"._ 








"3 


< 


-II 


* ^ 


-< 






* 2 














^ II 








" II 








H 






II .^ 








* £ 
£ 








■-J 






















K 






















V 
























.—1 


to 




t 
4^ 














9§ 


_ 


c^ 




DO 
















d 
o 


+J 




Is 














^ u 


■ H 


(U 




c^-^ 
















3 






Is 














a 


13 
d 
o 

U 








boa 
























1 



272 



a 
o 

a: 

I 






H 




K 




H 


^ 


K 


i^ 


O 




k. 


tb 


Q 


o 


^; 


-r^ 








O 


r/t 




O 


J^ 


M 


r« 


1" 


H 


m 


k 

t" 




H 


:ii 




< 


1^, 




O 

O 




a 


H 


'^ 


k. 


■^ 


n 


H 










H 


O 

H 


^ 






o 


o 


< 




hj 


n 


< 












M 






k 


Oh 


a 


w 


H 




ti 




O 




In 




T 




X 




O 












'^ 




w 




cS 




U-i 




>1 




W 

























a) 






-d 


a- 


lU 


m 


J=l 


^ 



aj_g 



o 



O =1 

« S 



o-( 






' 


J -J 










Tl 


OT! 


O 




a 






a 
"3 




73 




« 


■-; 


>- 


'- 




r; 


o 


-.-' 



- 


-r! 








S 


















































fl' 




O 


M 

^ 


•^ 


-n 


s 


U) 


Oj 



^■a 



Q^ ^ 






CO 




j= 


h(i 


Qj 


■t- 


-^ 








^ 


o 




a; 


S 


>. 








rt 




« 




-M 










d 




C) 




> 






o 


a 




a 






0) 

H 


N 


+3 


^ 


(J 




, 


tf 




H 


3. 



273 



274 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1955 

The qualitative behavior of the entries in Table V is easily understand- 
able. The Components /''„„o„*'' refer to the longitudinal magnetoconduc- 
tivity when both electric and magnetic fields are in one of the principal 
directions of the valley, the a direction. Since a magnetic held in such an 
a direction does not change the a component of the velocity of the car- 
rier, this longitudinal magnetoconductivity must vanish: 

I'\.J" = 0. 

It is easily verified that, for our model, the principal directions of a valley 
are the only directions in which the longitudinal magnetoconductivity 
contribution vanishes. Since the relative longitudinal magnetoconduc- 
tivity Act/ct is necessarily ^ 0, and for a cubic crystal is the negative of 
the relative longitudinal magnetoresistivity Ap/p, we can conclude that 
on our model the vanishing of the longitudinal magnetoresistance in any 
direction is possible, at least for cubic materials, only if the direction in 
question is a principal axis of all the valleys. It can further be shown, 
though we shall not give the details here, that lack of constancy of t 
over an energy shell, far from upsetting this conclusion, merely makes it 
impossible for the longitudinal magnetoresistance to vanish in ami direc- 
tion. 

The nonvanishing magnetoconductance effects can be described in 
terms of the current due to the force exerted by the magnetic held on the 
transverse Hall current. This current is proportional to the Hall current 
and inversely proportional to the effective mass — call it mj* — in the 
direction normal to the Hall current and to H, this being the direction of 
the force producing the second-order current. The Hall current, ajs we 
have noted in Section 6, is proportional to the zero-order current, hence 
to the reciprocal of the effective mass ms* in the direction of E, and to 
the reciprocal of the effective mass the- h* in the direction normal to E 
and H. 

To employ these ideas specifically, consider first the component 
o"BffjSfl*'', which measures the change in current in the a direction pro- 
duced by a magnetic field in the /3 direction. Here m^* = nia*, m^- n* 
= m^*{y 

(y ^ a, &) (72) 

the minus sign coming in because the second-order current is in the direc- 
tion opposite to E. When we insert the mas.s-dependence of o-nMH into 
(fi9) and combine with the /''^^p^''* of Table V, we do in fact find that 
CTaffSi?''' contains the masses only as indicated in (72). (The fact that 



a, 0), 7nj* = 


m.„*. Thus 




1 

cc ■ 


ma*~niy* 



THANSPOKT PROPERTIES OF A ilANY-VALLEY SEMICONDICTOU 27o 

Foasti''' tlcpencls on tfio masses in a much more complicated way is due 
to oin- choice of the dcHiiiiig ef|Uation for it, (69) : wc chose to write this 
piiuution so that it involved only the directly measurable ciuantities 
(Til iiiid nil , and the dimonsionlcss f|uantitics A and F^-^ag "■ This choice 
is tlie most con\'enicnt one for comparisons with experiment, but is less 
simple conceptually than a choice gi\'ing the factor (72).) The remaining 
independent component, F^ffas'''', may be analyzed similarly. It repre- 
sents the second-order current in the a direction due to an E in the ;3 
direction and an H in the direction midway between a and 0. Here 
!„^* = ms*, mn-H* = m^*(Y ^ a, ji), mj* = m„*. The Hall current is 
weaker by 2"' than for the previous case, because of the 45° angle 
between E and H, and since the force producing the second-order current 
is at 45° to the ^ direction, we must put a second 2^^''' in the denominator. 
Thus 



"^"^ 2ma*mg*my* 2mi*m3*m;i* 



(73) 



This, again, can be verified to follow from (69) and Table V. 

To apply these results to experimental magnetoresistance data in the 
region of proportionality to H', it is necessary, as has been mentioned 
above, to derive an experimental magnotoconductivity tensor from the 
observed magnetoresistance. For a cubic ciystal the magnetoresistivity 
tensor can be described by three constants h, c (not to be confused with 
velocity of light), d, defined by 

^^,fj^^^Jy^^,(jlMl±JlMl±JlE^ (74) 
pH' j- r 

where Ap is the change of resistivity p due to a small field H, and where 
the axes are those of the crystal. The equations relating the constants 
/>, c, d to the corresponding constants describuig the magnetoconductivity 
tensor have been given by Pearson and Suhl.^* From these equations the 
components of a^^^^ can be expressed in terms of the empirical constants 
h, c, d. The results are tabulated in the last row of Table VI. 

If the ratios of these components are compared with the ratios of the 
corresponding F^^«3 , one can check the correctness of an assumed model 
and determine the ratio m^f'/m^*. From the absolute values of the 
a^^,^^ one can then determine .1 . A further check is provided if data are 
available at more than one temperature, since the ma.ss ratio should come 
out roughly independent of T, while the variation of A with T should 

i«G. L. Pearson and H. Suhl, Phys. Eev., 83, p. 768, 1951. 
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accord with a reaaonable picture of the effects of impurity, inter-, and 
intra-valley scattering on the form of T-(Ae). 

An analysis of this sort has been carried out for n type silicon. For 
this substanc^e the longitudinal magnetoresistance nearly vanishes in 
directions of the type (100). From what has been said above, this almost 
requires that the band stnicture have valleys on the (100) axes in K- 
space, and that the rela.\ation time be practically a function of energy 
only. Fitting the remaining magnetoconductivity constants gives 
m*/mi^ ^ 5; this ratio comes out independent of temperature as it 
should. It agrees with the ratio determined by cyclotron resonance. 

9. HALL EFFECT AND MAGNETORESISTANCE FOR LARGE MAGNETIC FIELDS 

As the theory of Hall and magnetoresistance effects for large magnetic 
fields is rather complicated mathematically, it wiU suffice for our purposes 
merely to outline the approach which can be used and to quote a few 
results without proof. Some of the details can be found in the papers of 
Abeles and Meiboom^ and of Shibuya. 

To treat these we solve the transport equation (61) for the distribution 
function f'^' and calculate the current density from f'^'. This gives the 
electrical conductivity tensor (t^.(H), which can be inverted to give the 
resistivity tensor p^j, . The antisymmetric part of p^^ determines a Hall 
coefficient (in general slightly orientation -dependent)," and the sym- 
metrical part determines the magnetoresistance. 

The solution of (61) can be carried out either by summing the series 
(64), or directly by guessing that f'^' will be a linear function of the 
velocity components, with coefficients which are functions of energy. 
These coefficients can be determined by solving a set of three simultane- 
ous equations. 

As H -^ '-^ , the conductivity tensor ct^,(H) becomes singular, the con- 
tribution ff^/'' of the iih valley taking the form 

(o ne <Ae7-> HM, . , 

<y.. -^ TT- -TT^ ... .rr. , ...Irr. , ... *rT . ^ general 



Ny <A€> mi*Hi' + rm*H^^ -f- m^Hz^ 



o-o m^'^HM. 



Ny m^*Hi^ + m2*H-? + m^*Hz^ 



(75) 
for a cubic crystal 



where n is the total density of carriers, Ny the number of valleys, Ae the 
distance from the band edge, T(Ae) the relaxation time, o-q the conduc- 
tivity at H = 0, in'^^ the average inertial mass defined by (15), and 

" G. L. Pearson and C. Herring, Physica, to appear. 
" Dexter, Lax, Kip, and DresBelhaus, see Reference 5. 
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Hi , H2 , Hs are the components of H along the principal axes of the 
valley. Summing (75) on i gives a limiting a^y which is still proportional 
to H^Hy ; its determinant therefore vanishes, and it cannot be inverted 
to give the limiting p^„ . It turns out that (75) suffices to give the limiting 
value of the longitudinal magnetoresistance as H — > x , but that to get 
the limiting value of the transverse magnetoresistance it is necessary to 
evaluate contributions to o-^^(H) of order l/H^ The limiting Hall coeffi- 
cient can be obtained from the contributions of order l/H to (r;,„(H). 

The following points are noteworthy : (i) The limiting value of the Hall 
coefficient as H —> <xi is R = +(l/nec), independently of the arrange- 
ment and mass anisotropy of the valleys, and of the dependence of t on 
energy. Even if t is not a mere function of energy, the same limiting form 
obtains, (ii) For E parallel to H, the ratio cr{H — > «> )/ao depends on the 
arrangement and mass anisotropy of the valleys, but not on the law of 
variation of t with energy, as long as r is a function of energy only. 
(iii) For E not parallel to H, both (t{H -^ co) and a(H — * «)/o-o depend 
on the behavior of t(Ac) as well as on the band structure. 
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Appendix A 

THE RELAXATION TIME ASSUMPTION 

The use of the relaxation time concept represents a great simplification 
of the Boltzmann equation for all kinds of transport phenomena. To be 
strictly correct, one ought to describe the scattering processes which the 
carriers undergo by a transition probability S(K., K') defined as the prob- 
ability per unit time of a transition from an initial state K to a final state 
in unit volume of wave number space centered on K'. (Conservation of 
energy will usually limit these transitions to a certain surface in K'-space, 
so that S will contain a delta function of energy; however, this complica- 
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tion has no bearing on the remarks of this paragraph.) The Boltzmaiin 
equation, for the case of Maxwellian statistics, therefore takes the form 

^f{K) ^ (dKK)\ _^ f [/(K')S(K', K) - ,f(K).S(K, K')] dK'. (Al) 
dl \ dt /Fii-i.is ■' 

This is in general an integral equation for the distribution function /, or 
for the part E -f'^' of / which is linear in the electric field. Our task in this 
Appendix is to say a few words about the validity of approximating (Al) 
by an equation of the form 

dl \ dl /FU-uh \ r{K) / 

and to examine the validity of the further approximation t{K) = T(e). 
We shall give only a rather brief discus-sioii of these c|uestioiis, however, 
as a future publication" will give a more thorough treatment and include 
a discussion of the solution of (Al) when these approximations fail. 
To begin with, let us consider the special class of cases for which 

^(K', K) = S(K', K*) (A3) 

where K* is the state in the same valley as K but with opposite velocity. 
If (4) is inserted for / in the integral of (Al), the /"" term contributes 
nothing, and if (A3) is satisfied 

j'E-f'''(K').S'(K',K) r/K' = 0, (A4) 

since f"" is an odd function of velocity while S is even. Therefore, to the 
first order in LJ, (Al) reduces to (A2) with 

1/t(K) = fs(K,K') dK' (A5) 

For collision processes which do not satisfy the velocity-randomizing 
condition (A3) we may assess roughly the legitimacy of using a T(e) by 
comparing, for different initial states on the same constant energy sur- 
face, the mean rate with which scattering destroys the component of 
velocity in the original direction. ^Phis rate of loss of velocity defines a 
1/r which if isotropic is known to l)e usable in (A2) for energy-conserving 
scattering proces.ses in the simple moilel,''* and for the many-valley model 
it is a reasonable presumption, borne out l)y the more rigortms treatment 

"W. Shockley, Electrons and Holes in Semiconductors (Van Nostrand 1951) 

p, 251 et aeq. 
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of Reference 0, that if this rate of loss of velocity is nearly constant over 
ati energy surface, then (A2) is legitimate with r = T{e). 

We shall consider five types of scattering in turn, three by phonons 
and two l)y impurities. For the phonon processes, scattering of a carrier 
from K to K' involves absorption of a phonon of wave vector K' — K, 
or emission of one with the negative wave vector. The matrix element 
for such a process is of the form"" 

N '" l/r/Tr Tr'\ [absorption 

il/(K,K')^ _^^ ii^^^foA (Afi) 

(.V + 1) ' J [ emission 

where A'' is the occupation number of the lattice mode involved, w is its 
frequency, and f '(K, K') is proportional to the matrix element, between 
states K and K', of the perturbation of the electi'onic Hamiltonian due 
to a static displacement of the nuclei of the lattice, of unit amplitude in 
this mode. The scattering function .S(K, K') which wc have used above 
is given by 

.S(K, K') - y E I '1/{K. K') f 5U(K) - e(K') ± /ic] (A7) 

where S is the Uirac delta function and the summation is over absorption 
(upper sign) and emission (lower sign), and over the various branches of 
the vibrational spectrum. 

(i) Inkr-VaUcy Lattice Scattering. Here the magnitude of K' — K is 
large compared with the distance of either K or K' from the band edge 
point nearest it. Moreover, the total change in ^oj as K or K' ranges over 
a constant energy surface will Ije « hT if the energy I'elative to the band 
edge is ^'/;7\ Therefore the percentage variation of MCK., K') over such 
an energy surface will be .small. And unless the energy relative to the 
band edge is extremely small, the surface in K'-spacc which makes the 
argument of the delta function vanish will be nearly a constant energy 
surface. We conclude that »S'(K, K') can be taken to be independent of 
K and K' when either ranges over a constant energy surface in its valley, 
and in pai'ticulai', (A3) applies. Therefore inter-valley lattice scattering 
is described by a i-elaxation time which is given by (Ao) and is a function 
of energy only. 

(ii) Intra-Valley Scattering bij Optical Modes. For a nonpolar crystal 
this case is essentially the same as the preceding, since the matrix element 
C(K, K') is substantially etjual to its limiting value as K' — * K. For a 
polar material the longitudinal polar optical modes give a C(K, K') a 

^" See, for pxiim|)lp, Reference U), p. 520. 
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I K' _ K p\ In a cubic crystal this contribution is independent of the 
direction of (K' - K), but if the energy surfaces are anisotropic the 
variation of C(K, K') with | K' - K | will suffice to prevent 7-(K) from 
being constant over an energy surface. 

(iii) Intra-Valleij Scattering by Acoustical Modes. Here (A6) and (A7) 
again apply, but with the simplification that for all but extremely slow 
carriers the ftoj in tlie argument of the delta function can be neglected, 
since only very low-energy phonons are involved. However, C(K, K') 
need no longer be independent of direction. The reason for this is 
that according to the deformation potential concept,^^ C depends on the 
strain associated with the lattice mode in question. If only the volume 
dilatation affected C, as was the case for the sunple model, C would 
be independent of the direction of the phonon wave vector (K — K'}, 
since the dilatation amplitude in a compressional wave of unit dis- 
placement amplitude is independent of direction. But for a many- 
valley model both shears and dilatations can produce deformation 
potentials, and the nature of the shear stram in a shear wave of unit 
amplitude depends strongly on the direction of propagation. There- 
fore C may be a function of the direction of (K — K'), in any par- 
ticular valley. 

For a valley whose K''' lies on a threefold or fourfold symmetry axis 
of a crystal the deformation potentials due to the different possible types 
of strains can be expressed in terms of two constants, Zd , Su , which 
appear in the theory of piezoresistance, Eq. (CG) below, and which are 
defined thus: Let «r{r = 1 to 6) be the six components of the stram 
tensor, relative to the principal axes of the valley, the z axis bemg 
taken along the symmetry axis of the valley. Then a dilatation in the 
two directions normal to the symmetry axis (wi = 1/2 = m/2, Us = 0) 
produces a band edge shift Ed"- A uniaxial shear (id = U2 = — V2, 
U3 = u) produces a shift "EuU. As will be shown in detail in Reference 6, 
it is not hard to evaluate the dependence of C(K, K') on the direction 
of K' — K, in terms of Zd and Hn and the elastic constants of the crystal. 
Thus we can calculate the variation of S(K, K') over an energy surface 
by (A7), in terms of the constants Sj , E,. , and the anisotropy of the 
effective mass. 

We shall presently exhibit the results of some calculations, made by 
the procedure just outlined, of how much the rate of loss of initial 
velocity varies over an energy surface. However, we shall first present an 
argument that for given deformation potentials, the greatest variation of 
this quantity may be expected to occur for spherical energy surfaces, 

''^ See, for example, Reference 19, p. 520 et seq. 
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and that usually very little variation will occur for extremely prolate or 
oblate surfaces. In other words, for almost all values of Sd and Su , the 
use of a T(e) will be justified for very prolate or very oblate surfaces. Fig. 
11 shows the argument. In (a) is shown a spherical energy surface cen- 
tered on a band edge point on some symmetry axis in K-spaee. We 
expect the greatest difference in relaxation time r(K), defined in terms of 
rate of loss of initial velocity, to be that between a point Ki where the 
symmetry axis cuts the sphere and a point K2 90° around the sphere from 
Ki . Backward scattering processes for Ki , i.e., those which almost re- 
verse its velocity, have K' — Ki vectors almost parallel to the symmetry 
axis, while forward scattering processes, which take Ki to a state K' with 
almost the same velocity, have K' — Ki almost at right angles to the 
axis. For a carrier initially at Ka , on the other hand, the vector K' — Ka 
is almost perpendicular to the axis for backward scattering, while for 
forward scattering it may range from almost parallel (plane of paper) to 
perpendicular (normal to paper). For this case, therefore, a dependence 
of the scattering function on the inclination to the symmetry axis will be 
quite effective in producing an anisotropy of t(K). Now consider the 
situation for a very prolate energy surface, as shown in (b) of Fig. 11. 
This figure has been derived from that of (a) by a horizontal extension 
and a vertical contraction. Now all the dotted lines repre-sonting vectors 
K' — Ki or K' — Ka are nearly parallel to the symmetry axis, and the 
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Fig. 11 — Compiirison of intra-vallcy lattice scattering processes for a valley 
with spherical energy surfaces, (a) , and a valley with highly prolate surfaces, {b} . 
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corresponding scattering functions must therefore l)e very nearly the 
same. Of course, forward scattering from Ki through sufficiently small 
angles has K' — Kj nearly normal to the axis, and the same is true for 
backward scattering from K2 in a small range close to 180°. But it is 
clear that as the energy surface becomes more prolate these cases form a 
smaller and smaller fraction of the total of possible final states K'. Thus 
for extreme prolateness t(Ki) -* t(K2). For strongly ohhite energy sur- 
faces most of the K' — K vectors approach normality to the symmetry 
axis, and a similar conclusion holds. The argument fails if, and only if, 
the scattering probability practically vanishes for K' — K along the 
symmetry axis (prolate case) or perpendicular to it {oblate case). 

Fig. 12 shows the results of some calculations of t(K2)/t-(Ki) carried 
out by E. Vogt for the worst case, that of spherical energy surfaces, and 
for an actual case, that of valleys on a ( 11 1 ) axis with m,|*/mj.* = 1 .3/.08, 
the value found in cyclotron resonance experiments on n geiTnanium." 
The anisotropy of the elastic constants has been assumed to be that for 
germanium. For the spherical surfaces calculations were made for valleys 
centered on (100) and (111) axes, but the results were found to be in- 
distinguishable. A comparison of the full curve (spherical surfaces) with 
the dashed one (prolate surfaces) shows, as expected, that with a highly 
anisotropit^ effective mass the anisotropy of the relaxation time is much 
less than for the spherical case, except near the ratio S,//E„ = — 1 . This 
is the ratio for which modes with wave vectors along the symmetry axis 
arc incapable of scattering. We conclude that for intra-valley lattice 
scattering the assumption of a relaxation time dependent only on energy 
will fail over a considerable range of the deformation potential parame- 
ters if the effective mass is isotropic, but only over a moderate range 
near E,//Hu = — 1 if the effective mass is very anisotropic. 

(iv) ScaUrring by Ionized Impurities. To date the quantum theory of 
this effect has been developed only on the crude basis of treating the 
fluctuations of potential due to a random arrangement of ions as a small 
perturbation on the motion of the charge carriei-s.^^ The result is that the 
effective matrix element for scattering between two states K, K' on the 
same energy surface is a function of j K' — K | which is sharply peaked 
at very small values of this quantity, at least when the density of impuri- 
ties is not too high. This means that the principal effects come from small- 
angle scattering, a fact well known in the classical theory of impurity 



=2 H Brooks, Phya. llcv., 83, p. 8 & f), 1951 (eoiit.iiiiis U'poHi'iipliiciil error.s); 
C. Horic, Tflh. U. Sci. Rep. 34, 27 (1950); iiml C. HcrriiiK, uripuMiHhed. 
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Fig. 12 — Atiisolropy of the rehixiition time for intra-valley lattice scattering, 
as a function of the ratio of tlie deformation potential coefficient "Ed for two-dimen- 
sional dilatation In the coefficient Sn for uniaxial shear, as these quantities are 
defined in Appendix ,\ or in lu|uatioii {C6). The ordinate is the ratio, for the 
points K'j and K| of Fig. 11, of the effective relaxation time r defined by t--' — 
— (initial velocity)"' (time rate of change of mean velocity, <iue to scattering). 
The elastic anisolropy assumed is that for germanium. 



scattering.^*' Now the probability of a collision with a given range of 
I K' — K I i-s proportional to the sciuare of this matrix element and to 
the number of pow-sible final slates K' within the given range of distances 
from K and in a given .small range of energy. A little calculation .shows 
that this number is greater when K is near the Ki of Fig. 11(b) than when 
K is near K^ . Morcovei', the fractional loss of the velocity component in 
the initial direction Is greater for a collision at Ki than at K2 . The result 



^' E. CoiivvHI and V. F, Wrisskoi)f, Phys, Rev., 77, p. HS8, 1950; R. S. Cohen, 
L. Hpitzer, Jr., and V. M. Uoutiv, Phvs. T^ev,, 80, p, 230, 1950; and L. Hpitzcr, Jr. 
and R. Harm, Phvs. Rev., 89, p. 977, 1953. 
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is that carriers near Ki have a much shorter effective relaxation time for 
impurity scattering than do those near Kj . 

Thus if the energy surfaces are very anisotropic, the assumption t = 
T(e) may be expected to be a poor approximation when ionized impurity 
scattering is importaiU. 

(v) Neutral Impurity Scattering. In the simple theory the scattering of 
charge carriers from neutral impurities in hydrogen-lilce states is mathe- 
matically equivalent to the scattering of electrons from hydrogen atoms.^* 
At the temperatures at which such scattering is important the wave- 
length of the incident carrier is usually » the diameter of the wave 
function of the center, so s wave scattering predominates. Therefore the 
scattering is isotropic. For a many-valley model, however, the situation 
seems at first sight more complicated, since the effective mass is aniso- 
tropic, and the centers are not spherically symmetrical. However, it can 
be shown^ that, at least if the energy of the carrier is low enough, the 
scattered wave must be describable as an 5 wave in the space of the 
transformed coordinates defined by ^ 

^, = m,*'"x, (A8) 

In the corresponding momentum space (the ^-space of Appendix B) the 
surfaces of constant energy are spheres. It follows, that, at least at low 
energies, neutral impurity scattering satisfies (A3) and is describable by 
a constant relaxation time. 

Appendix B 
equality of the energy-shell averages 



- i^mi*Di^ Hm2*Vo^, }4m7*vi 

Choose coordinate axes along the principal axes of the energy surfaces 
in any given valley. For \ = 1, 2, 3, let 

<P. = v^{m,*r- = APx/(mx*)"^ (Bl) 

Then 

A6 ^ I e(AP) - 66 ! = M E Va' (B2) 

so that the energy surfaces are concentric spheres in ^J-space. The den- 
sity of ^-vectors consistent with the periodic boundaiy conditions is of 
course uniform, like that of AP. The average of J^^m^W over an energy 
shell is therefore the average of tp^ over a spherical shell, and thus ob- 
viously independent of X. 

" C. Efginsoy, Phya. Rev., 79, p. 1013, 1950. 
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Appendix C 
explicit calculation of elastoresistance constants 

Let <j^y be the conductivity tensor of a crystal, defined by the relation 

U - Z '^i^yE, (CI) 

between current density j and electric field E. Let Ua& be the strain tensor, 
defined by the relation 

5x„ = 2] Wa^x^ (C2) 



between displacement 6x and initial position x of any point in the body. 
Then the elastoresistance of any crystal is described mathematically by 
the fourth- rank tensor da^y/dUas , or more conveniently, for a cubic 
crystal, by the dimensionless tensor — o-~ da^y/dUag , where a is the scalar 
conductivity in the unstrained state (a^^ = aSf^v). The task of this appen- 
dix is to obtain this tensor by calculating the strain variation of the cur- 
rent contributions (7), which in terms of components take the form 



3^ 



Ci> 



= — E 

kT Ap(o,s 



rr(Ae'*')i;.,.X 



(C3) 



We shall base the calculation on assumptions (a) through (e) of Section 
5. Assumptions (b) and (d), regarding existence of a T(Ae) and nondegen- 
erate statistics, are already contained in (C3). According to assumption 
(a), we shall neglect any effect of strain on the relations between v, AP, 

and Ae''* = | e — e*'' |. Thus we may replace v^v, in (C3) by its average 



L^r over an energy shell, and treat this average for given Ae as unin- 
fluenced by strain : 



VuV, 



H^e'-' 



nil' 







—^ Q 










J-'UIP 



(C4) 



When the crystal is strained, the only things that we assume to change 
in (C3) are the population factor/*"' and the relaxation time T(Ae*'^). 
For some groups of electrons ^ i.e., values of i and e — the product 
t/*"' will be increased, and for others it will be decreased. If on the aver- 
age this product is increased in the valleys whose conductivity tensors 
are most favorably oriented to the direction of E, and decreased in the 
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others, the total current will be increased by the strain, and vice versa. 
When assumption (e) of Section 5 is fulfilled, so that two of the m*'s, 
say m* and m>*, are ef|ual, the tensor (C4) ta.kes the form 






mi* K'-'>-' Wa 



«3* nil*/ J 



(C5) 



L. 



(■) 



as is easily verified by inspection. Moreover, it is easily seen from sym- 
metry that the only strain components which can alter e''' in first order 
are the isotropic dilatation and the shear compounded out of an extension 
along K''* and a contraction in both directions at right angles. Mathe- 
matically expressed, we must have 

^ = EA. + ^.kJ^'K.'Vk''' (C6) 

where Sj and E„ are constants independent of the valley i, the subscripts 
referring to "dilatational" and "uniaxial" effects respectively. The elas- 
toresistance tensor is, for nondegenerate concentrations, 

1 dtj^r _ ^ v^ Y"* 

tn^ya^ = — - -—— — TTfT 2-, Zj 

(C7) 

(V) I q I 6 — Cf I _ yp ,.(0) "^ oe 

'" kT du.0 V de^'' dUas, 

where e^ is the Fermi level, T*'\e) is the relaxation time in the ith valley, 
and L;,/'' is given by (C5). The second term in brackets in (C7) repre- 
sents the effect of strain on the transition probabifity for inter-valley 
scattering. 

We shall now combine and simplify the e(|uations just given. The be- 
havior of the Fermi level is simple: by assumption (a) of Section 4 it 
does not shift in shear, and for extrinsic concentrations it shifts with the 
band edge in compression. Mathematically, 

^ - ^ E ^ - (h. -K ^y, (C8) 

dUag Nv i dUa0 \ 3 / 

where Nr is the number of valleys. By virtue of the fact that t^ is a func- 
tion only of Ae and the differences {e'*' — e*'''), 

and it is easily seen that when (CO) and (C8) are inserted into (C7) the 
Srf terms disappear. Since by symmetry dr'^/de"^ = d,('i)/de'\ (C9) 
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implies 

E 1^ = (CIO) 

Using this and (C6) and (C8) in (C7), and processing further by (C5), 
we get 



ynuvaa — 



e 



-J ii) jy- (J)ft^ (J)"l 



(Cll) 






1 _ 1 

wis* mi* 



where the upper sign is for n-type, the lower for p, n is the total carrier 
concentration, angular brackets with a subscript i or i, j mean averages 
on valleys i or i and j, and angular brackets without subscripts mean 
Maxwcllian averages as defined in Section 2. Substituting from (14) for 
th(^ conductivity a- = nen, we get finally 



niuycd ^-32,, 



/<:ita7"Vae'^ K/'^A^^" KJ''Ks'\ 



(012) 

{mi* - m.3*) 



^'' \ <AeT>~ A^'''- ■ K"'^^- /;.j1 {m,* + 2m;,*) 

As in (C7), the first term in stiuare brackets in (012) represents the 
effect of the strain ou the relative populations of the different valleys, 
and the second term represents the effect on the inter-valley scattering 
probabilities. We may note the following features: 

(1) The trace 2^„m^„ao \-aiiishes identically, by virtue of (C9) etc. 
This means that an isotropic dilatation produces no elastoresistance 
under the assumptions we are using. 

(2) The elastoresistance is proportional to the auisotropy of the effec- 
tive mass within a valley. 

(3) »ivi\i = mu vanishes for valleys of the (100) type, while mnn = 
niivi-i (niu = »(|o)for valleys of the (111) type. 

{'[) The elastoresistance is proportional to l^T at temperatures low 
enough for intei'-valley scattering to be fiozeii out. Moreover, when kT 
» the iiiter-valley hwiij) of Section 3, the variation of 6t, 'dtj with energy 
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Ac relative to the band edge is easily shown to be one of proportionality 
to t/Ac, if impurity scattering is negligible. Under these conditions 
(probably never achieved by Si and Ge in the extrinsic range) the elas- 
toresistance is again proportional to \/T, but with a larger factor of 
proportionaUty than at low T. 

NOTATIONS 

A dimensionless average involving relaxation 

time, Equation (68) 

B dimensionless function of mass anisotropy, 

Equation (39) 

C(K, K') factor in matrix element for a scattering pro- 

cess K — > K', Equation (A6) 

c velocity of light; magnetoresistance constant, 

Equation (74) 

Bij factor in matrix element for scattering from 

valley i to valley j, Equation (16) 

E electric field , 

e electronic charge 

F^^aQ factor determining the anisotropy of the mag- 

netoconductanco tensor Equations (70), 
(71) 

F^yoff'^^ contribution of the ith. valley to above, Equa- 

tion (69) 

/ distribution function for charge carriers. 

/'"' same in absence of perturbing fields, Equation 

(1) 
f'^' etc. change of / in perturbing fields. Equations (4), 

(31), (60) 
G dimensionless average involving relaxation 

time, Equation (57) 
H magnetic field 

k Planck's constant/27r 

j density of electric current 

K wave number vector for a charge carrier 

K''^ value of K for the ith band edge point (center 

of the I'th valley) 
k Boltzmann's constant 

L;,^''' average of Vf,Vy over an energy shell in the ith 

valley 
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M(K, K'), Mij matrix element for lattice scattering, Equa- 

tion (A6) or (16) 

m normal electron mass 

m\* effective mass in Xth principal direction of a 

valley 

m}'^ inertial average of the ?%*, defined by Equa- 

tion (15) 

m^rap elastoresistance tensor — o- da^JdUa^ 

]^^ JVa number of quanta (phonons) in a given lattice 

mode 

JVr number of valleys or band edge points 

n number of free charge carriers per unit volume 

P crystal momentum ftK 

Q^ electronic part of thermoelectric power 

q wave number vector for a lattice mode 

R Hall constant 

S(K, K') transition probability K to unit rfK' at K', 

Equation (Al) 

s spin quantum number of a charge carrier 

T absolute temperature 

i time 

u„3 , Ur strain tensor components 

V group velocity of a charge carrier 

Waiij, a) transition probability from valley i to valley 

j with absorption of a phonon of branch a 

Wciv, «) same but with emission 

Wi , W2 constants of dimension (time)~\ measuring 

coupling of carriers to intra- and intervalley 
modes, respectively, Equations (21)-(23) 

x^ y, 2 rectangular coordinates 

a (when not subscript) index labeling branches of vibrational spec- 
trum 

Y differential operator describing rotation of 

distribution by magnetic field, Equation 
(62) 

^g energy of a carrier relative to the band edge 

(Ae ^ 0) 

AP, AK departure of crystal momentum or wave vec- 

tor from valley center 

5,.„ ICronecker delta 
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S^ya antisymmetric coefficient — ±1 or 0, Equa- 

tion (33) 

be, 5v, etc. changes induced by strain 

£ energy of an electron 

tb band edge energy 

€*'* energy of the center of the ith valley (nor- 

mally = et) 

O Debye temperature 

K dielectric constant 

Ij. drift mobility 

fi,i Hall mobility 

Ed , S„ deformation potentials, Equation (C6) 

{ transformed coordinate, equation (A15) 

p resistivity 

(7 conductivity 

ffn conductivity in the absence of magnetic fields 

or other perturbations 

<T^.a , (T;,yas coefficients in the expansion of conductivity in 

powers of H, Equation (65) 

iT^ya/' contribution of I'th valley to magnetoconduc- 

tivity 

T relaxation time of charge carriers 

r''' ditto for carriers in valley i 

<p transformed wave vector, Equation (Bl) 

CO angular frequency of a lattice mode or an rf 

field 



